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Abstract 

One presents here the *^-Bernouhi- Taylor* formula of a new sort with 
the rest term of the Cauchy type recently derived by the author in 
the case of the so called V'-difference calculus which constitutes the 
representative for the purpose case of extended umbral calculus. The 
central importance of such a type formulas is beyond any doubt - 
and recent publications do confirm this historically established experi- 
ence. Its links via umbrality to combinatorics are known at least since 
Rota and MuUin-Rota source papers then up to recently extended by 
many authors to be indicated in the sequel. 

KEY WORDS: umbral calculus, Bernoulli formula, Graves-Heisenberg- 
Weyl algebra(**) 

AMS S.C. (2000) 05A40, 81S05, 01A45, 01A50, 01A61 

* see below : a historical remark based on Academician N.Y.Sonin arti- 
cle published in Petersburg in 19-th century. We owe this information and 
the article to Professor O.V.Viscov from Moscow. 



** see: C. Graves, On the principles which regulate the interchange of 
symbols in certain symbolic equations, Proc. Royal Irish Academy 6(1853- 
1857), 144-152 
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1 One Historical Remark 



Here are the famous examples of expansion 



oo 



,n-l 




n\ dx' 



n=l 



or 



oo 



n! dx ' 



^o = E(-i) 



n 



n=0 



where do is the divided difference operator while eo is at the zero point 
evaluation functional. If one compares these with "series universalissima" 
of J.Bernoulli from Acta Erudicorum (1694) (see commentaries in H2|) and 
with 



then confrontation with B.Taylor's "Methodus incrementorum directa et in- 
versa" (1715), London; entitles one to call the expansion formulas considered 
in this note "Bernoulli - Taylor formulas" or (for n oo) "Bernoulli - Tay- 
lor series" [1]. 

Information: Johann Bernoulli was elected a fellow of the academy 
of St Petersburg. Johann Bernoulli - the Discoverer of Series Universalis- 
sima was "Archimedes of his age" and this is indeed inscribed on his 
tombstone. 

2 Introduction 

A. From here now ^/^ denotes an extension of (^)n>o sequence to quite 
arbitrary one (the so called "admissible"- see: Markowsky) and the specific 
choices are for example : Fibonomialy -extended {{Fn) - Fibonacci sequence ) 
(■f3-|)n>o or just "the usual" (^)n>o or Gauss g-extended {:;^)n>o admissible 
sequences of extended umbral operator calculus - see more in [1611131 [HI [T5] . 

The simplicity of the first steps to be done while identifying general 
properties of such [6,7,8,16,13-15] -(/'-extensions consists in natural notation 
i.e. here - in writing objects of these extensions in mnemonic convenient 
upside down notation [15] , [H] 




dx ' 



fc=0 



V'(n-l) 



! = n^{n - 1)^1, n > 0,x^ 



^p{x — 1) 
tp{x) 



(1) 
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x| = x^{x - l)^{x - 2)^...{x -k + l)^ 



(2) 



x^{x - l)^...{x -k + l)^ 



^(x — l)il){x — 2)...il){x — k) 
ip{x)ilj{x — l)...^p{x — k + 1) 



If one writes the above in the form x^ 




one sees that the name upside down notation is legitimate. 

You may consult [10] , [9] ,[14,15] for further development and use of 
this notation. 

With such an extension we frequently though not always may " ^-mnemonic' 
repeat with exactly the same simplicity and beauty most of what was done 
by Rota \10\ [T5] . Accordingly the extension of notions and formulas with 
its elementary essential content and context to general case of ip- umbral in- 
stead of umbral or (/-umbral calculi case only - is sometimes automatic [lOj , 

[T5] (see corresponding earlier references there and necessary definitions) . 

B. While deriving the Bernoulli- Taylor ■(/'-formula one is tempted to 
adapt the ingenious Viskov's method [2] of arriving to formulas of such type 
for various pairs of operations. In our case these would be -(/^-differentiation 
and ^-integration (see: Appendix). However straightforward application of 
Viscov methods in ■(/'-extensions of umbral calculus leads to sequences which 
are not normal (Ward) hence a new invention is needed. This expected and 
verified here invention is the new specific *^ product of analytic functions or 
formal series. This note is based on [3] where the derivation of this new form 
of Bernoulli- Taylor *^ - formula was delivered due to the use of a specific 
*^ product of formal series. 

3 Classical Bernoulli-Taylor formulas with the rest 
term of the Cauchy type by Viskov method 

Let us consider the obvious identity 



in which ([4]) we now put ak = a^h^]a,h € A. A \s an associative algebra 
with unity over the field F=R,C. Then we get 



n 




(4) 



fc=0 



n 



J2aHl-ab)b'' = l-a"+^&^ 



,n+l. 



a,b e A 



(5) 



fc=0 
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Numerous choices oi a,b £ A result in many important specifications of ^ 
Example 1. Let denotes tlie linear space of sufficiently smooth func- 
tions f : F — > F. Let 

'f{t)dt, 



a: T — 




(a/)(x) 


h:T - 




(V)(^) 


l:T - 







\^f){^)\ (6) 
ax 



Then [b,a]=l-ab=eQ where is evaluation functional on T i.e. 

e^{f) = /(a) (7) 

Using now the text-book integral Cauchy formula {k > 0) 

(a'^f)(x) = l^^^-^^^f{t)dt, (8) 

and under the choice ([H) one gets from ^ the well-known Bernoulli- Taylor 
formula 

/(^) = E + (9) 

k=0 

with the rest term Rn+i{x) in the Cauchy form 

Rn+l{x)= ^ ^fin+l)it)dt (10) 

Example 2. [1] Let J-" denotes the linear space of functions / : Zj^ — > 
F;Z+ = NU {0}. Let 



x-l 
A:=0 

f(x + l)-f(x), (111 



It is easy to see that [b,a]=l-ab=eo where eo is evaluation functional i.e. 
eo(/) = /(O). 6 = A is the standard difference operator with its left inverse 



a : 


■■z+- 




{af){x) 


b ; 




— ^; 


ibf)ix) 


I : 






mix) 
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definite summation operator a. The corresponding A - calculus Cauchy 
formula is also known (see formula (31 p. 310 in [5j); 

(a'Dix) = ^^7— T^l=/W; k > (12) 

r=0 

where x— = x{x — l){x — 2)...{x — n + 1). 

Under the choice (jlip one gets from ([5]) the A - calculus Bernoulli - Taylor 

fomuia yy 

k=0 

with the rest term Rn-\-i{x) in the Cauchy A form 

Rn+iix) = Y: ^""'T ^" (A"+V)(r); (14) 

r=0 



4 " *^ realization" of Bernoulli identity. 

Now a specifically new form of the Bernoulli- Taylor formula with the rest 
term of the Cauchy type as well as Bernoulli- Taylor series is to be supplied 
in the case of -^-difference umbral calculus (see [5-8] and [9,10] and references 
therein). For that to do we use natural -^-umbral representation [13,14] of 
Graves-Heisenberg-Weyl (GHW) algebra [11,12] generators p and q and then 
we use Bernoulli identity (fTSl) 



n 



ki n! 



fc=0 

derived by Viskov from ^ under the substitution (see (28) in [2]) 

ao = 0, afc = (-l)'=(g)'=- V(fe - 1)!, A: = 1, 2, ... 
due to pq^ = q^p + nq'^^^ (n=l,2,...) resulting by induction from 

[p,q] = l (16) 



Example 1. The choice p = D = and q = x — y,y G F;xf{x) = xf{x) 
after substitution into Bernoulli identity (jlSp and integration dt gives the 
Bernoulli - Taylor formula Q. 
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Example 2. The choice [2] p = A and q = xoE^^ where E"' f{x) = f{x+a) 
after substitution into Bernoulh identity (|15p and "A - integration" X]r=o 
gives the Bernoulh - Mac laurin formula of the following form (a, a; E Z, y = 
1 — i?"^) with the rest term Rn+i{x) 

n ^ 

= E If (-l)'^'(vV)(a) + i?n+i(a); (17) 

fc=0 

r=0 

Example 3. Here f^''^ = d^f and f{x)*^g{x) = f{x^)g{x) - see Appendix. 
The choice p = and q = {z = x — y) where x^x" = j^^y^x^~^^ after 
substitution into Bernoulli identity (fTSl) and "0^ - integration" J^d^t (see: 
Appendix) gives another Bernoulli - Taylor ■i/'-formula of the form: 

n 
fc=0 

with the rest term Rn+i{x) in the Cauchy-form 

^n+l (X) = r dqt{x - tr*^ (t)(it (20) 

In the above notation x'^*^ = 1, x"^*^ = x *^ (2;("-i)*</') = x *^ ... *^ x = 
^x"; n > 0. 

Naturally d^x"^*^ = nx^'^~^^*^ and in general f,g- may be formal series for 
which 

d^{f H S) = {D f) *^ g + f *^ {d^g) (21) 
i.e. Leibniz *^ rule holds [T3 | [14 } [TS]. 

Summary: These another forms of both the Bernoulli -Taylor formula with 
the rest term of the Cauchy type [3j as well as Bernoulli - Taylor series are 
quite easily handy due to the technique developed in [13l [I^ where one may 
find more on *^ product devised perfectly suitable for the Ward's "calculus 
of sequences" p] or more exactly *^ is devised perfectly suitable for the so- 
called V' - extension on Finite Operator Calculus of Rota (see [9l [lOl [HI [15] 
and references therein) 



6 



5 Appendix 

*^ product 

Let n — ijj = ijjn', ipn ^' n> 0. Let be a linear operator acting on formal 
series and defined accordingly by d^x'^ = n^x"'~^. 

We introduce now a intuition appealing S^-difference-ization rules for a spe- 
cific new *^ product of functions or formal series. This *^ product is what 
we call: the ^-multiplication of functions or formal series as specified below. 
Notation A.l. 

X *^ = x^{x'^) = -^^^-x^^^'i n > Q hence a; *^ 1 = {l^)~^x ^ x 

therefore x *^ a\ = al *^ x = x a = a x = a{l^)~^x and Vx, a G F; 
f{x) *^ x^ = f{x^)x'^. 

For k x"*^ x*^ 7^ x^*^ x"' as well as x"*^ x*^ 7^ x"+'^ - in general. 

In order to facilitate the formulation of observations accounted for on the 

basis of ^-calculus representation of GHW algebra we shall use what follows. 

Definition A.l. With Notation A.l. adopted define the *^ powers of 
X according to x"** = x ^^x^"^'^^** = x^{x^^~^^*'t') = x *^ x *^ ... *^ x 
= ^x"; n > 0. Note that x"** *w, x''** = J^x^'^+''>^ ^ x^** *w, x"** = 
^ {n+k)*^ for /c 7^ ra and x^** = I. 

This noncommutative ^-product *^ is devised so as to ensure the fol- 
lowing observations. 

Observation A.l 

a) 5^x"** = nx("-i)*'/'; n > 

h) exp^[ax] =exp {ax^}l 

c) exp [ax] *^ (exp^{/3x^}l) = (exp^{[a /3]x^})l 

d) d^{x^ *^ x'^*^) = {Dx^) *^ x"** + x^ *^ {d^x^*'!') 

e) d,^{f *^ g) = (Df) *^ g + f *^ {d^g) ; f,g - formal series 

f) f{x^)9{x^) 1 = f{x) *v six) ; g{x) = g{x^)l. 

Umbral "~" Note: g{x) = g'(x^)l defines the map ^. g ^ g i.e. 
P ^ P which is an umbral operator. 

T/j-Integration Let: c^ox" = x""^. The linear operator do is identical 
with divided difference operator. Let Qf{x)f{qx). Recall also that to the 



7 



"dq difference-ization" there corresponds the (/-integration which is a right 
inverse operation to " g-difference-ization" . Namely 



F (z) ■■^(U) (z) := (1 - g) z f; (q'z) q 



(22) 



I.e. 



F (z) ^(^l^^yz) = {l-q)z (^f; q'^Q'^^ {z) 



Of course 



as 



{'^-<l)z-^^]{z). (23) 
l-qQ 



dgo = id (24) 



q 



^ '^%({l-q)z^-^]=id. (25) 



(1-9) ' l-qQ. 
Naturally (j25p might serve to define a right inverse operation to " g-difference- 
ization" (dgip) (x) = jf^doip{x) and consequently the " g- integration" as 
represented by ()22p and (|23p . As it is well known the definite g-integral is 
an numerical approximation of the definite integral obtained in the g — > 1 
limit. 

Finally we introduce the analogous representation for difference-ization 

= n^do] n^x"""^ = n^x^~^; n > 1 (26) 

Then 

[ x''= (x—] x'^ = - — 3-— n > (27) 
and of course = / 

d^o [ =id (28) 

Naturally 



d^o I f{t)d^t = f{x) 
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The formula of "per partes" -i/^-integration is easily obtainable from (Obser- 
vation A.l e) and it reads: 

pb j'b 

/ {f*i'di.gmd^t=[{f*^g){t)t- {Df*^g){t)d^t (29) 



Closing Remarks: 

I. All these above may be quite easily extended [E] to the case of any Q G 
End{P) linear operator that reduces by one the degree of each polynomial 
|16j . Namely one introduces |15j : 
Definition A. 2. 

XQ G End{P),XQ : F[x] ^ F[x] 

such that (x") = ^^+1) ^ 0; where Qqn = nqn-i- Then ★g product 

of formal series and Q-integration are defined almost mnemonic analogously. 

II. In 1937 Jean Delsarte [17] had derived the general Bernoulli- Taylor 
formula for a class of linear operators 5 including linear operators that reduce 
by one the degree of each polynomial. The rest term of the Cauchy-like 
type in his Taylor formula (I) is given in terms of the unique solution of a 
first order partial differential equation in two real variables. This first order 
partial differential equation is determined by the choice of the linear operator 
8 and the function f under expansion. In our Bernoulli- Taylor -formula (16)- 
(17) or in its straightforward -kQ product of formal series and Q-integration 
generalization - there is no need to solve any partial differential equation. 

III. In [18] (1941) Professor J. F. Steffensen - the Master of polynomials 
application to actuarial problems 



(see : http://www.math.ku.dk/arkivet/jfsteff/stfarkiv.htm) 



supplied a remarcable derivation of another Bernoulli- Taylor formula with 
the rest of "Q-Cauchy type" in the example presenting the "Abel poweroids" 

IV. The recent paper [19] (2003) by Mourad E. H. Ismail, Denis Stan- 
ton may serve as a kind of indication for pursuing further investigation. 
There the authors have established two new q-analogues of a Taylor se- 
ries expansion for polynomials using special Askey- Wilson polynomial bases. 
As "byproducts" their important paper includes also new summation theo- 
rems, quadratic transformations for q-series and new results on a q-exponential 
function. 

V. Let us also draw an attention to two more different publication on the 
subject which are the ones referred to as [20, 21]. The (/-Bernoulli theorems 
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are named here and there above as g-Taylor theorems. The corresponding 
{q, /i)-Bcrnoulli theorem for the 9g /^-difference calculus of Hahn [22] might 
be also obtained as the the one (q, 0)-Bernoulli i.e. (/-Bernoulli theorem 
constituting here the special case the Viskov method [2] application. This 
is so because the Sq^^-difference calculus of Hahn [22] may be reduced to 
Q'-calculus of Thomae- Jackson [5,23] due to the following observation . 
Let 

he R, {Eq^hLp){x) = ip{qx + h) 

and let 

(P> ,n\(^\ - ipjqx + h) 

= (30) 

Then (see Hann [5] and [22]) 

dq,h = ^dqE-. _h^. (31) 

'1-q '1-9 

Due to (30) it is easy now to derive corresponding formulas including Bernoulli- 
Taylor 5g_/j-formula obtained in [24] by the Viskov method [2] which for 

9^ 1,/i^O 

recovers the content of one of the examples in [2] , while for 

1,/i^ 1 

one recovers the content of the another example in [2] . The case /i — >^ is 
included in the formulas of g-calculus of Thomae- Jackson easy to be specified 
: see [22] (see also thousands of up-date references there). For Bernoulli- 
Taylor Formula (presented during PTM - Convention Lodz - 2002) : contact 
[23] for its recent version. 

The comparison of the all above quoted ways to arrive at extended 
Bernoulli formulae we leave for another exhibition of similar investigations. 

VI. As indicated right after Observation A.l.e) the rule g{x) = g{x^)l 
defines the map g ^ g which is an umbral operator P ^ P. It is 
mnemonic extension of the corresponding q - definition by Kirschenhofer 
[24] and Cigler [25]. This umbral operator (without reference to to [24,25]) 
had been already used in theoretical physics aiming at Quantum Mechanics 
on the lattice [26]. The similar aim is represented by [27,28] (see further ref- 
erences there) where incidence algebras are being prepared for that purpose 
(Dirac notation included). As it is well known - the classical umbral [29,30] 
and extended [10] finite operator calculi may be formulated in the reduced 
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incidence algebra language. Hence both applications of related tools to the 
same goal are expected to meet at the arena of GHW algebra description of 
both [14,13]. 
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